A geometrical approach to the dynamics of spinor condensates I: Hydrodynamics 
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In this work we derive the equations of motion governing the hydrodynamics of spin-_F spinor 
condensates. We pursue a description based on standard physical variables (total density and su- 
perfluid velocity), alongside 2F 'spin-nodes': unit vectors that describe the spin F state, and also 
exhibit the point-group symmetry of a spinor condensate's mean-field ground state. The hydrody- 
namic equations of motion consist of a mass continuity equation, 2F Landau-Lifshitz equations for 
the spin-nodes, and a modified Euler equation. In particular, we provide a generalization of the 
Mermin-Ho relation to spin one, and find an analytic solution for the skyrmion texture in the incom- 
pressible regime of a spin-half condensate. These results exhibit a beautiful geometrical structure 
that underlies the dynamics of spinor condensates. 



I. INTRODUCTION 

A central theme of contemporary atomic physics exper- 
iment is the dynamics of Bose-Einstein condensates and 
other correlated atomic gases. Of particular interest are 
mixtures of several species such as Fermi-Bose mixtures 

S 111 J and bosonic gases with an internal spin degrees 
of freedom, i.e., spinor condensates. Spin-one and spin- 
two spinor condensates have been realized as particular 
hyperfine states of alkali atoms 0, H, § . In addition the 
trapping and cooling of ^^Cr atoms has led to the real- 
ization of a spin-three condensate Q • On the theoretical 
front, since the initial work of Ohmi and Machida Q and 
Ho [9j numerous interesting works have followed that dis- 
cuss ground states, dynamics, and topological excitations 
of such systems (see, for instance, a review in [loj). 

Recently, it was shown within mean-field theory that 
the ground states of spinor condensates exhibit a high 
degree of symmetry. This symmetry is opaque in the 
standard spinor description of the condensate. On the 
other hand, the symmetry is transparent in the so-called 
reciprocal state representation. Here, one uses the fact 
that the mean field ground state of a spin-_F condensate 
can be described by 2F coherent spin states orthogonal 
to it. Each one of these so-called reciprocal states is 
fully spin-polarized, pointing along some direction on the 
unit sphere. Since there are 2F such reciprocal states, 
the ground state is uniquely described (up to an overall 
phase) in terms of the 2F points on the unit sphere [ll| . 
For typical spinor condensate Hamiltonians, these points, 
which we denote "spin nodes", form highly symmetric 
configurations. For instance, an = 2 condensate has 
a cyclic phase, where the spin nodes are arranged in a 
tetrahedron, as well as a square phase. 

The spin-node description of the ground states of 
spinor condensates provides an intuitive geometrical de- 
scription of the state of the condensate. In addition, 
it provides a parametrization which readily exhibits the 
hidden point-group symmetries of the state. Despite its 
appeal, however, this parametrization has not been used 
to describe the dynamics of spinor condensates. Our goal 



in this paper is to provide a complete description of the 
hydrodynamics of spinor condensates in terms of such 
spin-nodes. 

The description of the dynamics we provide in this 
paper is hydrodynamic in the sense that it focuses on 
the low energy dynamics of the system associated with 
locally conserved quantities, or with the slow elastic de- 
formation of spontaneously broken degrees of freedom. 
As with any fluid, such a macroscopic hydrodynamic de- 
scription is highly illuminating and useful. Here we de- 
rive such a description using the density, superfluid ve- 
locity, and the spin-nodes (the 2F vectors on the unit 
sphere) as our basic degrees of freedom. In addition to 
the Euler equations, which describe mass, momentum, 
and energy conservation, we obtain 2F Landau-Lifshitz 
equations for the dynamics of the spin-nodes. Further- 
more, our derivation gives a natural generalization of the 
Mermin-Ho relation which connects the vorticity in a fer- 
romagnetic spinor condensate with the Pontryagin index 
of the order parameter. 

The paper is organized as follows. Sec. [Til provides 
general background on spinor condensates, and reviews 
recent progress on the hydrodynamic description of fer- 
romagnetic condensates [12[- In Sec. IIIII we present the 
spin node representation of spinor-condensate degrees of 
freedom, and derive several useful identities within this 
formalism. In Sec. IIVI we proceed to obtain hydrody- 
namic equations for the spin-half condensate, using the 
spin-node formalism, and find an analytic solution for a 
skyrmion configuration. In Secs.lVland lVIl we derive the 
general hydrodynamic equations of motion for the spin- 
one condensate, and then for an arbitrary spin-i^ conden- 
sate, which includes a generalization of the Mermin-Ho 
relation p^ . 

The treatment of the spin degrees of freedom in this 
paper is exact, and it accounts for the full geometrical 
structure of the hydrodynamics of spinor condensates. 
But the precision of the hydrodynamic description here 
comes at a price: this formalism becomes increasingly 
complex as the spin F grows, and, for large F, the anal- 
ysis of its exact form becomes impractical. Nevertheless, 
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the equations derived here even for large F become quite 
useful in their linearized form. In an accompanying pa- 
per we show how approximate methods - such as 
linearizing the equations of motion about mean-field so- 
lutions - elucidate the low energy properties of spinor 
condensates with arbitrary spin in a powerful and ele- 
gant way. 

II. BACKGROUND 
A. Hydrodynamics of spinless BECs 

For a single component BEC, it is natural to expect 
a simple hydrodynamic description in terms of density 
and flow velocity. We take the time-dependent Gross- 
Pitaevskii equation (GPE) as our starting point: 

idt^ = -^W^i^ + gpij (1) 

where ip = y^e** is the macroscopic wave function, and 
p = j-i/jp is the density (here and after for notational 
simplicity we will use scaled units). This equation can 
be recast into a the form of local momentum and mass 
conservation laws; with the superfluid velocity v = V0, 
one obtains [isj : 

dtp = -V(pv) ; Av = -V (^gp - (2) 

where I?( = (9f + v • V is the material derivative. The 
first of these is the mass continuity equation, while the 
second is the Euler equation for a fluid, where a quantum 
pressure term appears. 



Lamacraft's approach assumed an incompressible liquid 
with a wavefunction restricted to the ferromagnetic phase 
(assuming large g and C2) 

V-a = e*'^$a(n) (5) 

where $a(n) is the highest eigenstatc of n • F. A sub- 
stitution of this wavefunction into Eq. yielded the 
following set of hydrodynamic equations |12l | : 

V • V = (6) 

An = -n X V^n. 

2 

Once the density is eliminated, we notice that the spin 
dynamics are given by a Landau-Lifshitz equation with 
the material derivative Dt = dt + v ■ . In addition, the 
vorticity is related to the Pontryagin density by 

F 

V X V = —EaPjUa (Vn/3 X Vti-^) . (7) 

This identity is widely-known as the Mermin-Ho rela- 
tion d, [l^. Among other things, such a relation has im- 
portant consequences for the topological defects in ferro- 
magnetic condensates d [13, [2ij2i . Such hydrodynamic 
equations were also derived in [30| to describe magnetic 
properties of quantum hall systems. 

Making use of this simple description, Lamacraft 
showed that the helical configuration of the ferromag- 
netic condensate [l^ is unstable. In general, it is clear 
that such a geometric description simplifies, at least con- 
ceptually, the analysis of spinor-condensate dynamics. 

C. General magnetic ground state of spinor 
condensates, and the reciprocal state representation 



B. The hydrodynamics of ferromagnetic BECs and 
the Mermin-Ho relation 

In a series of recent experiments, the quench dynam- 
ics of a ferromagnetic spin-one condensate was explored 
[l6l . [TtI . ITSi . .19] . These experiments motivated Lamacraft 
to develop a hydrodynamic framework for the ferromag- 
netic BEC in terms of the superfluid velocity, and the di- 
rector of its ferromagnetic order n [l^l . This description 
is particularly illuminating when considering the instabil- 
ities of the system. Th is p roblem was also theoretically 
considered in Refs. 0, M, M, HI, [H, [H, [2^ . 

The GP Lagrangian density describing such a ferro- 
magnetic spinor condensate is given by 

C = ii'ldti^a - ^V^: • VV'a ' \gp^ - \c2P^ni^ (3) 

where a = — i^, . . . , -F is summed over all F^ eigenstates, 
and 

P-E^aV-a ra=-Y,ra^a,A- (4) 

a ^ ah 



A general spin-_F spinor-condensate is described by a 
macroscopic wave-function with 2F -\- 1 complex compo- 
nents ipa (a = —F, . . . ,F). When quantum fluctuations 
are unimportant, the condensate dynamics is described 
by the time-dependent Gross-Pitaevskii equation 

^9ti'a=-lvH'a + ^^ (8) 

where V is the spin-dependent interaction energy. The 
interaction energy is given by the set of parameters gs, 
with 5 = 0, 2, ... , 2F describing the two-particle interac- 
tion strength in the S total angular momentum channel: 

^nt = lY.9srarb {ab\Sm) {Sm\a'b') i^a'^b'- (9) 

In the above, {ab\Sm) are Clebsch-Gordan coefficients. 
Note that this expression can also be written as the ex- 
pectation value of an operator: 

Vint = l ,{ij\Smnt\i'\\i^\ (10) 



where 

Vint = as ISm) {Sm\ = ^ gsVs- (11) 

S,m S 

In this expression, Vs projects into the total spin S scat- 
tering channel. 

The classical (mean-field) ground states occur for uni- 
form condensates which minimize Eq. ^ for fixed den- 
sity p. This minimization was carried out for F = 1 
F = 2 [m, and F = 3 [13, [H yielding a multitude of 
magnetic phases, which minimize V\at in different regions 
of {gs} parameter space, only one of which (for every F) 
is ferromagnetic. 

Indeed, quite generally, a spin-F spinor condensate 
may exhibit several flavors of paramagnetic rather than 
ferromagnetic behavior in its ground state. For example, 
a spin-one condensate may exhibit the so-called nematic 
phase, where ipi = ip-i = 0, and i/'o ~ 1- The expecta- 
tion value of the magnetization for such a state is clearly 
zero along any direction, (F • n) = 0. But in the ab- 
sence of a ferromagnetic director, n, can we still describe 
a spinor condensate's magnetic state geometrically? 

Such a geometrical method was put forward in Ref. 
based on the use of spin-coherent states. A spin- 
coherent state |$n) is the eigenvector of the opera- 
tor F • n with the largest eigenvalue. The method of 
Refs. [ll| relies on finding the set of 2F spin-coherent 
states, { I which annihilate the ground state of 
a uniform condensate: 

($n.K'Gs)=0. (12) 

The 2F states |$ni) provide (up to an overall phase) a 
unique description of the magnetic spin-state of the con- 
densate at each point in space. Such reciprocal spinors 
give a natural generalization of the ferromagnetic direc- 
tor to the case of paramagnetic condensates. Instead of 
the geometrically opaque 2F + 1 complex numbers tpa , it 
allows a description of the magnetic state in terms of 2F 
unit vectors, n-i, or points on the unit sphere. 

In addition to its geometrical transparency, such a de- 
scription also reveals the highly symmetric nature of the 
mean-field ground states. All the paramagnetic phases 
found so far correspond to a spin node configuration 
which is invariant under point symmetry group opera- 
tions, and sometimes under a larger symmetry. The ne- 
matic phase of the F = 1 condensate, for instance, is 
described by two antipodal spin nodes. F ~ 2 conden- 
sates can exhibit a nematic phase as well, but also a phase 
in which the spin-nodes are the vertices of a square, and 
a phase with the spin nodes at the vertices of a tetrahe- 
dron. Such phases are illustrated in Fig. [T] 

The reciprocal-spinor description was so far only uti- 
lized to discuss equilibrium properties of spinor conden- 
sates. The remarkable geometrical properties and hidden 
symmetries of the mean-field ground state, however, pro- 
vide ample motivation for employing the spin nodes to 
obtain a complete description of the dynamics of spinor 
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FIG. 1: Possible phases for the spin-one and spin- two con- 
densates. Spin-one condensates have ferromagnetic (la) and 
nematic phases (lb) while spin- two condensates have ferro- 
magnetic (2a), uniaxial nematic (2b), square biaxial nematic 
(2c), and a tetrahedral (2d) phases. 

condensates. In the following sections we will develop the 
tools necessary for such a description, and use them to 
derive both a hydrodynamic description, as well as small 
oscillation dynamics near mean-field ground states. 

III. SPIN NODE DESCRIPTION OF SPIN-F 
MAGNETIC STATES 

The reciprocal-spinor states so far define the spinor 
condensate's state only implicitly through Eq. In 
order to be able to use these variables directly, we must 
invert the relationship, expressing the spinor condensate 
Lagrangian directly in terms of these variables. In order 
to find this direct representation, we first separate the 
wave function into a piece corresponding to the overall 
density and phase, and a piece describing the local spin 
state. We write 

■Ipa = IpXa (13) 

where Xa is a normalized spin-F spinor 

Ex:xa = (xix) = 1 (14) 

a 

and the supcrfluid density is 

p=\^f. (15) 

A. Symmetrized spin-node representation 

As discussed above, a spin-F spinor \x) can be de- 
scribed by 2F reciprocal states. On the other hand, such 
a state can also be described by a fully symmetrized col- 
lection of 2F spin-half states. Each spin-half state can 
be parametrized in terms of two coordinates f2 = (9, (j)) 
on the unit sphere, 

\n) = cos (^) e^<^/2 IT) + sin (^-) e^^/^ |i) . (16) 
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In this representation, 
1 /^^ 



\i=i 




(17) 



where A/" is a normahzation constant, and the sum over a 
runs over the (2f[)! permutations of the 2F labels for the 
spin-half parts [SJI- In Eq- (HZl) we also defined as 
the (unnormalized) sum over permutations of the tensor 
product. 

The properties of the above formulation are most easily 
understood using the Schwingcr Boson construction |35| 
(for review, see [36[). Schwingcr Bosons provide an easy 
way to construct the Hilbcrt space of a spin-i^ spinor 
state. We define two Schwingcr boson creation operators: 
(x\w . An a) boson adds 1/2 to both the total spin, and 
to Fz; whereas a 6^ boson adds 1/2 to the total spin, but 
lowers Fz by half. In this notation 



Ftot — 


-(ata- 








i(atSH 






Fy = 






(18) 




-{a^a- 


- Pb). 





A spin-half spinor is written as: 

\n)=u\^)+v\i) = {ua^ 



vb^)\0) 



with |0) the Schwinger-boson vacuum. Here, u and v can 
be written in terms of the coordinates on the unit sphere 
as w = cos(6l/2)e''>/2 and v = sm{e/2)e-'">'/^ . 

A symmetrized tensor product of 2F spins within the 
SB formalism is simply written as: 



\n) = \ni...n2F) = 



2F 
i=l 



rbA |0) 



(19) 



with Ui and Vi parametrized in terms of 9i , as shown 
above. We refer to this collection of the 2F spin-half 
states which construct \ as "spin nodes". 

If we wish to calculate wavefunction overlaps using the 
Schwingcr Boson formalism, we can use Wick's theorem 
to obtain 



t2F 



(oin 
E n (oi { 

2F 
{a} »=1 



(a)* ~ 

u- a 



a + V. 



2F 



u - 'a' 



(b)-t 
u)j/a' 



-vf^b^ 
bA |0) 



|0) 



ib)\ 



(20) 

where tr is a permutation of the 2F indices that mark 
the spin-half parts. This result could have also been ob- 
tained directly from Eq. P7|) . Nevertheless, we find it 
instructive to demonstrate the simple Schwingcr Boson 
construction to obtain the symmetrized states. 



B. Connection between spin nodes and reciprocal 
spinors 

Since the symmetrized spin-node representation can 
be used to express any spin state directly, it makes a 
grossly overcomplete basis. Nevertheless, its usefulness 
arises since it perfectly reflects the spin-nodes formalism 
of the spinor-condensates ground states [111, [13] . In the 
following we will introduce the necessary new notation for 
the spin-node formalism; we summarize the new notation 
in Appendix El 

It is simple to see that a spin-coherent state can be 
written in terms of Schwingcr boson states as 



= (ua^+vb^j |0) 

.n) . 



(21) 



\n. 



Thus a coherent state can be thought of as 2F spin-nodes 
pointing in the same direction. (For a summary of the 
notation see Appendix 1X1) 

As described in Sec. Ill CI a reciprocal spinor is a 
coherent state [(fir)^^) = \il.r---^r) orthogonal to a 
given spinor |r2). Using the construction in terms of 
symmetrized spin nodes, we can write an equation to 
determine the reciprocal spinors for a particular state 

\n} = \ni...n2F): 



2F 



{{nrf^\n) = (2F)! Yl (i^rl^i) = 0. 



(22) 



This equation has 2F solutions, each corresponding to a 
different term in the product vanishing. That is, the ith 
solution of Eq. ^ is 



where 



(23) 



(24) 



is the time- reversed spinor of with 6^ ~ tt — 6i, <f>l ~ 
+ TT. Here we use the fact that a spin-half spinor is 
orthogonal to its time- reversed counterpart. 

The set of reciprocal coherent states is, therefore, just 
the set of states antipodal to the directions in the spin- 
node set {ft,i}i£i on the unit sphere. Using this construc- 
tion, Eq. provides us with a direct expression for a 
spinor condensate state in terms of its reciprocal vectors. 



C. Time derivatives of spinors and spin- nodes 

As stated above, our goal in this paper is to extract 
the dynamics in terms of individual spin nodes In 
order to do so, we must be able to isolate the dynamics of 
each spin node within |r2). Consider the time derivative 
of jJl) which will appear in the CPE. We can express it 
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as a sum of terms in which the time derivative operates 
on individual spin nodes: 

dt \n) = \dtn1n2 ...) + ...) + .... (25) 

The trick that aUows us to isolate individual spin nodes 
consists of taking the inner product of dt |r2) with the 
ith reciprocal-state of which is [(fi*)^-'^). All terms 
which do not involve a time derivative of identically 
vanish and we are left with the single term 

2F 

{{n^^^\dt\n) = {2F)i{nl\dtn,) [] {nl\n,). {26) 

We will make extensive use of this method for isolating 
the dynamics of individual spin nodes in the following 
sections. 

D. Geometrical parametrization of the spin- half 
components: Moving from 1^2;) to n^ 

All results above were concerned with breaking a spin- 
F spinor into its 2F spin-half parts, and with the 

correspondence between these spin-half parts and the re- 
ciprocal coherent states. We would like, however, to un- 
derstand the dynamics of spinor condensates not only in 
terms of the spinors, but also in terms of the unit 
vectors that describe them, n^, where is highest value 
eigenvector of F • . 

The first step in finding the equations of motion in 
terms of the spin-directors n^, is to establish an orthonor- 
mal triad (ej;,ej,,n) that parameterizes the space on 
in the vicinity of n^. In the following we will only con- 
sider a single spin-half part, and therefore we drop the 
index i. 

The first element of the triad is n itself: 

n = 2{n\F\n) (27) 

where F is the spin operator, acting in the spin-half 
Hilbert space. To complete the triad, we again use the 
time-reversed spin-half ket, jfi*) where (ri*|ri) = 0. With 
this, we can construct states pointing in the "x" and "y" 
directions with respect to n as 

\n,) = ^{\n) + \n')) (28) 
\ny)^^^{\n)+z\n')). (29) 

These states allow us to complete the orthonormal triad 
by defining 

e, = 2{n,\F\n,), ey^2{ny\F\ny). (30) 

From these we can construct 

e± = Gj; ± iey. (31) 



It is useful to note that F • e± act as raising and lowering 
operators. That is, 

F -6+ |f7) = ; F • e+ |f7*) = , (32) 

with similar relations holding for lowering operators. 

Note that there is an ambiguity in such coordinate sys- 
tems since e^; and can together be rotated about n 
which corresponds to the gauge choice for the spinors. 
That is, the gauge of a spinor can be changed by jfi) 
e'^ |ri) without changing its direction n. In general, 
quantities which are gauge invariant cannot depend on 
the parameterization of the spin and will only involve the 
unit vectors rij. We will adhere to the following gauge 
convention for parameterizing spin-half spinors: 

\n) - cos(0/2)e'^/2 It) + sin(0/2)e-'^/2 |i) (33) 

In this gauge-choice it is easy to sec that: 

ex = ey = LP (34) 

where 9 and (p are unit vectors from the spherical coor- 
dinate system. 

To complete the discussion, we make two observations 
that will simplify the following analysis. First, we express 
F in the basis of our triad as 

F = (F • n)n + (F • e^^e^ + (F • ey)ey 

= (F.n)n+i(F.e+)e_ + i(F.e_)e+ (35) 

In addition we note that we can use the spin operator 
F as a projection onto \ and its time-reversed partner 
m by 

\n){n\ = l + n-F (36) 

and 

|f2*)(r2*| = i-n-F. (37) 

E. Derivatives of spin-half spinors in terms of the 
triad (e^,, e^, n) 

The relations derived and recalled in the previous sec- 
tions allow us to also write derivatives of spinors in terms 
of vector quantities and their derivatives. The terms that 
we will encounter arise from terms such as the isolated 
time derivatives in Eq. (|26p . Let us now find this decom- 
position in terms of the triad (e^,, e^, n) and its differen- 
tial forms. 

Our goal is thus to find: 

a„ = i (17|(9„r2) and (f7*|a„fi) (38) 

in terms of (e^;, ej,,n) and their derivatives. We define 
aa in this form for reasons that will become clear later. 
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The first object in Eq. p8|) can be found by considering 
the quantity 

da{{n'\e_-F\n))^dal = 0. (39) 

Ahowing the derivative to operate on the bra, the ket, 
and the vector e_, we find 

(aarj*|fl*) + (n|9„rj) = -ie+ • 9„e_ (40) 

On the left-hand side we used the facts that F • e_ |fi) = 
and F e = On the right-hand side we 

used the fact that (ri*| F = ^e^ which can be verified 
from Eq. (|35|1 . It is easy to verify that 

from which we find 

which is the desired result. 

To obtain {^^\doS^) we use a similar trick. Starting 
with 

= (fl*|rj*) {n^\VL) = (17*1 i - n • F (42) 

we find 

a„((r!*|i-n.F|17)) = 0. (43) 

As before, allowing the differentiation to act on the bra, 
the ket, and n results in 

{n*\dM) ^ {^^\ {daxi) ■F\n) (44) 

where the term with (Safi*! vanishes since (^— n-F) jfi) = 
0. Now, using again the decomposition in Eq. (|35p . we 
readily find 

{n'\d,M) = \e+-d,,n. (45) 

This concludes all the tools we will need for our analy- 
sis below. We have found how to directly write a spin-F 
spinor in terms of its spin-nodes, and extract terms hav- 
ing to do with individual spin nodes out of sums arising, 
e.g., from differentiation. Furthermore we translated the 
spin-half representation of to a set of 2F triad bases 
(eij,, e^j^, Hi) which will allow us to parametrize the spin 
state geometrically. Appendix ?? summarizes the various 
notation introduced throughout this section. 

IV. HYDRODYNAMICS OF SPIN-HALF 
CONDENSATES 

One of our main goals is to write the exact (mean-field) 
equations of motion for a spinor condensate in terms of 



the spin nodes and the superfluid velocity and density. In 
this section we achieve this goal for spin-half condensates. 
The equations of motion can be trivially generalized to 
general spin-F condensates restricted to the ferromag- 
netic state, when the spinor \x) is restricted to be a co- 
herent spin-state. In this case the equations of motion 
for the condensate reduce to those we find below. 



A. Gross-Pitaevskii Lagrangian 

In this section we consider the Gross-Pitaevskii La- 
grangian. We begin by writing the Lagrangian in a re- 
vealing form, using the representation of the bosonic field 
which separates the spinor order parameter into a prod- 
uct of a density piece and a spin piece, i}}a — V'Xa- The 
GP Lagrangian is then: 

£ = ilpldtiJa - • '^V'a - Vint- 

= irdti^ + pat ~ i|(-iV - a)7/.|2 - ipT - Mnt- 

(46) 

where VSnt is the spin- related interaction and p = 
with -0 a complex field. Eq. defines the spin vector 
potential: 

at = I {x\dtx) ; a = i (x|Vx) , (47) 
and the quantity 

T = {d^x\dc.x) - (x|5aX) (9aXlx) • (48) 

An interesting observation is that the quantity T for a 
general spin F = N/2 can be identified with the CP^ 
model from quantum field theory [ssj . Notice that there 
is a U{\) gauge freedom in the density-spin decomposi- 
tion: 

^ ^ e^V, |X> ^ e-'^ |X> • 

The quantity T, however, is gauge independent. We 
make a gauge choice when we write the normalized |x) 
as in Eq. with written as Eq. ([TH). This 

U{\) gauge freedom is also reflected in an ambiguity in 
the choice of the triad arising from the spin-half parts of 
Ix) oc jn), since for each spin-part e^, and can together 
be rotated about n. The choice of a particular triad is 
set by the gauge choice. In general, quantities which are 
gauge invariant cannot depend on the paramctrization of 
the spin, and will only involve the unit vectors n^. The 
vector potential can be related to the superfluid velocity 
by 

v=4^(CVVa-V'aVC)=V0-a (49) 

pZl 

where 6 is the argument of t/i. So far we have not used 
the fact that the spin \s F ~ 1/2. 
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B. Geometric representation of hydrodynamic 
quantities 

Now that we know the quantities of interest in the 
spinor description of the GP Lagrangian, we can trans- 
late them to the hydrodynamic variables of density and 
magnetization direction. 

The most important quantity appearing above is the 
vector potential as defined in Eq. (|T7l) . Following the 
discussion in Sec. IIIIDI we see that for a spin-half con- 
densate the vector potential is 



condensate is 



(50) 



The analogy between a and the vector potential appear- 
ing in the Maxwell equations compels us to consider the 
antisymmetric field tensor f^^ = daap — dpaa- Through 
a series of manipulations this can be written purely in 
terms of n 



1 



(51) 



-{daey ■ n){dfjej: ■ n) - -(S/je^ • n){daex ■ n) 

^{ey ■ dan){ea: ■ dpn) - i(ey • di3n){e.j, ■ dan) 

1 
2 



-{ey X e^;) • {dan x dpn) 



= --n • {dan X dpn). 

Note that in the above we have repeatedly used the fact 
that V • 5v = for any unit vector v. The result is 
the Pontryagin topological density, which is the object of 
the celebrated Mcrmin-Ho relation for spin-half spinors 

[il,[3i]. 

The only remaining term is the gauge invariant quan- 
tity T, defined in Eq. For a spin-half state, we 
find: 



T = {dan\dan) - {dan\n) {n\dan) 

1 



(52) 



\{dan) ■ {dan). 



(e„ • dan){e+ ■ dan) 



Thus the T term signifies the stiffness of the superfluid 
with respect to magnetic gradients (as opposed to simply 
U{1) phase gradients). Also, since we identified T with 
the Lagrangian density of a CP^ model, we now reaffirm 
its equivalence with the nonlinear sigma model (ssj . 



C. Equations of motion for spin-half condensates 

Now that we clarified how the hydrodynamic variables 
arise in the GP Lagrangian density, we are ready to ap- 
proach the GP equations of motion. In terms of the orig- 
inal variables, the time-dependent GPE for a spin half 



idfipa = -^V^Vv + gp'4'a 



(53) 



where we note that the interaction energy for this case is 

1 



^nt = ■7:9P 



(54) 



Following the substitution tjja = ^'Xa, with Xa the entries 
of the spin-half spinor and contraction with {x\ we 
find 

idtip + i'at = ^(-«V - a)2-0 + ^T^j + gpip 

where = {at , a) is the vector potential introduced pre- 
viously. Substituting tp = /e*^ and multiplying both 
sides of the equation by e~*^ gives 

idtl -dtOf + fat - i(-V2/-z/V-v-2ivV/ + /«2) 



1 



rf + 9Pf 



(55) 



where v = V0 — a. The imaginary part of this gives 

dtp ^ -V • (pv) (56) 

which is a mass conservation equation. On the other 
hand, taking the real part gives 



ivy 1, 

2 / 2 



dtO + ^v^-at^^^-i-T-gp. (57) 



We take the gradient of both sides of this equation (using 
the identity V{v'^) = 2(v • V)v 2v x (V x v)) to get 



Av = e+(vx b)-V(.gp+iT-^^ 



(58) 



In this we have defined the "electric" and "magnetic" 
fields e and b in the usual way from the vector potential. 
That is, Cq, = f^t and bo, = (V x a)^ = ^ea/j^f/s^, with 
a, (3 and 7 indicating space directions, and the f tensor 
defined below. Also, note that we have used the material 
derivative Dt = 9* -I- v • V. The "electromagnetic force" 
appearing in the right-hand-side of the Euler equation 
is a new feature that is not present in single component 
condensates. This new type of quantum pressure arises 
from non-uniform spin textures in spinor condensates. 

Now we move on to find the equations describing the 
spin dynamics. To do this, wc contract the GPE with 
the time reversed spinor (xt|- This gives 

1 Vf 

« {xt\Dtx) = -2(2*a • (xtlVx) + 2-^ ■ (xt|Vx) 

+ (XtlV'x))). (59) 

Using the spin identities developed in Sec. Illll the follow- 
ing relations can be derived (xt|V^x) — ^da{e+ ■ dan) 
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and {dax\x) {Xt\daX) = ~\dae+ ■ dan.. In terms of vec- 
tors, the above equation is then 



An 



1 /l 



2 



V^n 



daf 



dan 



which can be rewritten as 
1 



pDtn = - (n X daipdan)) 



(60) 



(61) 



which is a Landau-Lifshitz equation. 

Thus, collecting everything, we can write down a com- 
plete set of equations describing the dynamics of a spin 
half condensate: 

~dtp V • (pv) 
V X V = -b 

Av = e + (v X b) - V (^gp + ^{danf - 

pDtn = i (n X daipdau)) 

where e and b are related to the spin direction through 
the field tensor 



fa/3 



/ 
















bz -by 




H 


-bz 


b. 


v 




by 


-bx 



-n • {daU X dfsn). 



(62) 

It is interesting to compare these results with those 
obtained previously in the incompressible regime p^ . 
Eqns. ([6]). Wc find that lifting the incompressibility con- 
straint leads to an Euler equation with effective electric 
and magnetic fields given by the Mermin-Ho relation. In 
addition, the superfluid density now enters the Landau- 
Lifshitz equation. 



D. Application: skyrmion texture 

As an example of the efficiency of the above hydro- 
dynamic equations of motion, let us consider a specific 
calculation: skyrmion textures in ferromagnetic conden- 
sates. For a standard U{1) vortex, the superfluid ve- 
locity close to the vortex core diverges as 1/r. For a 
scalar condensate, this causes the superfluid density to 
be depleted in a small region of order of the coherence 
length around the core. This can be energetically costly 
if the condensate is near the incompressible regime. On 
the other hand, this situation can be circumvented for a 
spinor condensate. Consider for example, a two compo- 
nent (spin-half) condensate (V'f ,'!/'i)i f-i^d take the J, com- 
ponent to have a U{1) vortex. Then around the vortex 
core, the density of ipi can be transferred to the vortex- 
free ipi keeping the total density across the vortex core 
finite. This is known as the skyrmion configuration which 
has been argued to be the relevant topological defect for 
ferromagnetic condensates [2^ . 




FIG. 2: A skyrmion configuration corresponding to Eq. (j66 



Let us now derive the analytic time-independent so- 
lution of the equations of motion in the incompressible 
regime having the skyrmion texture shown in Fig. [S] To 
this end, wc take the incompressible limit [T^ of the equa- 
tions of motion for the spin-half condensate obtained in 
Sec. IIV CI Neglecting z-dependence, these are: 



Ox Vy dy Vx 



V • V = 
1 



n • {dxn X dyii) 



An = -(n X V^n). 



(63) 
(64) 

(65) 



With small modifications, these equations can also be 
shown to describe the dynamics of condensates confined 
to the ferromagnetic phase of arbitrary spin in the incom- 
pressible regime. Our aim is to find stationary solutions 
of these equations having a skyrmion texture given by 

im 

n = (sin(/3) cos{(p), sin(/3) sin((p), cos(/3)) (66) 

where ip is the azimuthal angle and /3 is a function of r 
which is subject to the boundary conditions /3(r = 0) = 
and f3(r — R) — tt where i? is a distance far from the 
skyrmion center. Such a spin configuration is shown in 

Fig.m 

Given the form n in Eq. Eqns. (|63l64p can be 

solved to obtain the velocity profile. One finds 



sin^(/3/2) „ 
V = ip. 



(67) 



Note that the boundary condition /3(0) = suppresses 
the velocity at the origin which diverges as 1/r for the 
standard U{1) vortex. With the assumption of a static 
configuration, Eq. (|65|) reduces to 



(68) 



V • Vn = -(n X V^n). 
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With the expression for v in Eq. ((67|) . Eq. ((681) leads to 
the following second order differential equation for (3 



(69) 



With the boundary conditions, the solution of this dif- 
ferential equation is 



(Or) = 4tan~^ {r/R) . 



(70) 



This expression, along with the velocity in Eq. ([57)1 and 
the spin direction in Eq. (|66p constitute an analytic 
stationary solution to the equations of motion for the 
skyrmion configuration. 



V. HYDRODYNAMICS OF SPIN-ONE 
CONDENSATES 

A. Geometrical representation of spin one 
hydrodynamic quantities 

Now we move on to considering the more complicated 
case of the spin-one condensate. The spin-one spinor can 
be broken down into its two spin-half components and 
be written as |x) = /^/{fl\n) where = jf^il^a), 
where we again make use of the large-spin notation de- 
fined in Eq. The normalization factor for this case 
is found to be 



describes their coupling. This expression was previously 
obtained in Refs. [4l| and [i^l , where a geometrical rela- 
tion for the Berry phase of a spin-one spinor was given. 
The field tensor corresponding to this vector potential 
can also be similarly computed. The most simplified form 
we find is 



fa/3 



(3 + ni •n2)2 
(2ni • (9c.ni x dfjUi) + 2n2 • {dan2 x 9/3n2) 
+ (ni + n2) • (Saiii X di3n2 + dan2 x dptii)) 



(75) 



This is a generalization of the Mermin-Ho relation to 
the spin-one case. To our knowledge such an expression 
has not been previously derived. While its geometrical 
interpretation is not as immediate as the spin-half case 
(which is the Pontryagin density), this expression might 
be of use in computing topological invariant quantities 
for spin-one fields. This formula has a simplified form 
when locally restricted to mean-field ground states. For 
instance for the ferromagnetic sate (n = rii = n2) the 
above expression reduces to 



faf3 = -n • {daVi X d/^n). 



(76) 



It is also useful to note that for the nematic state (n = 
ni — — n2) the field tensor identically vanishes, fap = 0. 

Finally, the gauge invariant quantity T can be worked 
out to be 



= — I — ni • n2. 
2 2 



(71) 



It is also instructive to calculate the spin operator's 
expectation value. To start we can expand into products 
of spin-half expectation values 

(n|F|fl) = (rii|F|f)i) + (r!2|F|r!2) (72) 

+ {ni\n2) {n2\F + (1)2^1) f 1^2) . 

Then using the identity in Eq. (j36p . and the fact that 
in the factored expression, F is only acting on spin-half 
states, we obtain F jfJ) = ni -|- n2. Dividing this by 
the normalization, we get the spin-one expectation value 
of the magnetization: 



m= (x|F|x) =2 



ni -I- n2 
3 -f ni • n2 ' 



(73) 



By similar techniques, the vector potential for the spin- 
one case, with some work, can be written as 



1 



1 



aa = « {x\dx) = 2^iy ' (^a^^ix + -jf^2y ■ dae2x (74) 
^ 1 (n2 X m) ■ dgiii + (ni x n2) ■ dan2 
2 3 + ni • n2 

One sees that the first two terms in this expression are 
the vector potentials from the individual spin-half com- 
ponents while the final term, which is gauge invariant, 



T = 



r(5Qni • 9„ni + 9„n2 • dan2 + daiii ■ 80,112 



(3 + ni • n2)2 
ni • n2 da-ni ■ dan2 ~ ni • dc^2 n2 • i9ani). (77) 



This is an explicit representation of the CP^ model which 
can be viewed as a generalization of the nonlinear sigma 
model. Here, too, it is instructive to consider what this 
expression reduces to when locally restricting to mean- 
field ground states. For the ferromagnetic state, one finds 



T = -9an • 9o,n. 
2 



(78) 



On the other hand, for the nematic state T reduces to 



T = -a„n • a„n. (79) 



B. Spin-one condensate equations of motion 

We now proceed to do a similar analysis for the spin 
one problem. For this we note that the spin one GP 
energy functional has the form 



(80) 



where m is the expectation value of the spin-one op- 
erator. The first two hydrodynamic equations - the 
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mass continuity equation and the modified Eulcr equa- 
tion - are obtained, as before, by contracting the Gross 
Pitaevskii equation with (x| ■ The analysis proceeds along 
similar lines as the spin half case. However, for this case 
we need the generalization of the Mermin-Ho relation for 
spin one given in Eq. ()75p to give the field tensor and 
thus the effective electric and magnetic fields, in addi- 
tion to the spin one expressions for T Eq. (|77p and the 
magnetization m, in Eq. ([73]) . With these quantities, the 
first two equations of motion are 



dtp = -V • (pv) 



(81) 



and 



Av = e + (v X b) - V (^gp + C2pm^ + - 

(82) 

Next, let us discuss the spin dynamical equations. 
These are obtained by contracting the GPE with 
and (r2|r2||. As before, this causes several terms to van- 
ish since these spinors are orthogonal to |x) . Contracting 
with gives the following equation which gives the 

time derivative of the first node 



iei+-Ani = -irJ,^ei+-aQni-iei+-V^ni+C2/5ei+-m. 

. (^^^ 

A similar equation for the time derivative of n.2 is ob- 
tained by contracting with (il2^2l- In the above, we 
have collected the following terms into the parameter 



= 2 



y/p 3 + Hi • rij 

Wi ■ daUj - iUi ■ (rij X daUj) 



1 - n; • Uj 



+ i 



(rij X rii) • daiii + (rij x n^) • danj 
3 -t- rii • Uj 



(84) 



(85) 



two equations of motion, the mass continuity equation 
and the Euler equation are found, as before, to be 



Finally, separating the real and imaginary parts as F^ = 
(F^)' -I- i(F^)", we obtain the Landau-Lifshitz equations 

(A+^(Fi2)"a„)ni = inix((Fi2)'a„ni+V2ni)-C2pnixm 

(86) 

(A+^(r2i)"a„)n2 = in2x((F2i)'a„n2+V2n2)-C2pn2xm 

(87) 

This provides a complete set of equations describing the 
dynamics of the spin-one condensate. 



VI. HYDRODYNAMICS FOR GENERAL 
SPIN-i^ CONDENSATES. 

Now that we have considered the hydrodynamic equa- 
tions for spin-half and spin-one condensates in detail, in 
the following we will consider the general case. The first 



dtp = -V • (pv) 



and 



Dtv = e -f (v X b) - V 



2Kn 



-T - 



2VP 



(88) 



(89) 



The effective electric and magnetic fields again fol- 
low from the field tensor fajs constructed from Oq = 
i(x\9aX)- For a general spin, however, such quantities 
are cumbersome to express directly in terms of the spin 
nodes, and we will refrain from doing so. 

To obtain the Landau-Lifshitz equations, we contract 



the GPE with {{nl) 



i{nl\dtn. 



n2F 



Doing this gives 



-da log 



(90) 



In this expression, we have used the notation for interac- 
tion energy introduced in Eq. (jlOp . In addition we have 
introduced the quantities A,;, 



A, 



= i2Fy.l[{n,\nl) . 



(91) 



Instead of writing Eq. [501 in terms of vectors as in the 
previous sections, we will stop at this point. This equa- 
tion provides a natural starting point in the analysis of 
the linearized equations of motion which will be devel- 
oped in the companion paper Ref . [l3| . 



VII. CONCLUSIONS 

One of the most striking and surprising features of 
spinor condensates is the hidden symmetry of their mean 
field ground states. In this work, we have strived to bring 
this symmetry to the forefront, and describe hydrody- 
namics of these condensates using the objects that ex- 
hibit the hidden symmetry, namely, the spin-nodes, and 
reciprocal coherent states. 

In this work, wc derived the hydrodynamic equations 
of motion for condensates of general spin, demonstrated 
their use in the computation of the skyrmion configura- 
tion of a ferromagnetic spin-half gas, and generalized the 
Mermin-Ho relation to spin-one condensates. 

A hydrodynamic description tries to capture the low- 
energy behavior of a continuous medium in terms if con- 
served and other simple quantities. Therefore, it comes 



11 



particularly handy when we consider spinor-condensates 
close to their mean-field ground state. In a companion 
paper [3] , we will concentrate on small oscillations of the 
spinor condensate in the vicinity of the mean field-ground 
state. As we will show, it is there that the hidden point- 
group symmetry becomes most apparent and accessible. 
Using the spin-node formalism, and the parametrization 
of the spin-nodes in terms of a stereographic projection, 
we will reduce the problem of finding the 2F spin-wave 
eigenmodes to a simple question of decomposing a repre- 
sentation of the appropriate point symmetry group to its 
irreducible representations. We will also provide a sim- 
ple recipe that allows the direct extraction of the con- 
densate's spin-wave eigenmodes using the derivatives of 
the spherical harmonics, coupled with the knowledge of 
atomic orbital degeneracy lifting under a crystal field. 



Acknowledgments 

It is a pleasure to acknowledge useful conversations 
with E. Demler, T.-L. Ho, I. Klich, A. Lamacraft, and 
especially A. Turner. We would like to acknowledge 
the hospitality of the KITP, supported by NSF PHY05- 
51164. We are also grateful for support from the Sherman 
Fairchild Foundation (RB); the Packard and Sloan Foun- 
dations, and the Institute for Quantum Information un- 
der NSF grants PHY-0456720 and PHY-0803371 (GR); 
and CIFAR, NSERC, and CRC (DP). 



APPENDIX A: NOTATION 

In this Appendix, for convenience, we collect in one 
place the notation used for the representation of the 



spinors. Normalized spinors of arbitrary spin arc denoted 
by Ix). Non- normalized symmetric combinations of spin- 
half states are denoted as (with bold fonts) 



\n) = \nin2...n2F) 



(Al) 



This is then related to |x) by 



\x) 



(A2) 



Coherent spin states occur when all of the spin-half con- 
stituent spins point in the same direction. We denote 
these by 



2F 



(A3) 



Wc next define the spin state corresponding to |0) with 
its ith component time-reversed. We denote these by 



|T,n) = |f7ii72...a*...i^2F) 



Finally, we define the projection operator V to be 



(A4) 



v = i-\x) (xl 



(A5) 



[1] F. Schreck, L. Khaykovich, K. L. Corwin, G. Ferrari, [9 

T. Bourdel, J. Cubizolles, and C. Salomon, Phys. Rev. [10 

Lett. 87, 080403 (2001). 
[2] Z. Hadzibabic, C. A. Stan, K. Dieckmann, S. Gupta, [11 

M. W. Zwierlein, A. Gorlitz, and W. Ketterle, Phys. Rev. 

Lett. 88, 160401 (2002). [12 
[3] G. Roati, F. Riboli, G. Modugno, and M. Inguscio, Phys. [13 

Rev. Lett. 89, 150403 (2002). [14 
[4] J. Stenger, S. Inouye, D. M. Stamper-Kurn, H. J. Mies- [15 

ner, A. P. Chikkatur, and W. Ketterle, Nature 396, 345 

(1998). [16 
[5] H. Schmaljohann, M. Erfiard, J. Kronjager, M. Kottke, 

S. van Staa, L. Cacciapuoti, J. J. Arlt, K. Bongs, and [17 

K. Sengstock, Phys. Rev. Lett. 92, 040402 (2004). 
[6] M. S. Chang, C. D. Hamley, M. D. Barrett, J. A. Saner, [18 

K. M. Fortier, W. Zhang, L. You, and M. S. Chapman, 

Phys. Rev. Lett. 92, 140403 (2004). [19 
[7] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and 

T. Pfau, Phys. Rev. Lett. 94, 160401 (2005). [20 
[8] T. Ohmi and K. Machida, J. Phys. Soc. Japan 67, 1822 

(1998). ' [21 



T. L. Ho, Phys. Rev. Lett. 81, 742 (1998). 
M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, 
A. S. De, and U. Sen, Advances in Physics 56, 243 (2007). 
R. Barnett, A. Turner, and E. Demler, Phys. Rev. Lett. 
97, 180412 (2006). 

A. Lamacraft, Phys. Rev. A 77, 063622 (2008). 

N. Mermin and T.-L. Ho, Phys. Rev. Lett. 36, 594 (1976). 

R. Barnett, D. Podolsky, and G. Refael, arxiv:0903.1863. 

C. J. Pethick and H. Smith, Bose-Einstein Condensation 
in Dilute Gases (Cambridge University Press, 2002). 

L. E. Sadler, J. M. Higbie, S. R. Leslie, M. Vengalattore, 
and D. M. Stamper-Kurn, Nature 443, 312 (2006). 
S. R. Leslie, J. Guzman, M. Vengalattore, J. D. Sau, 
M. L. Cohen, and D. M. Stamper-Kurn, arxiv:0806.1553. 
M. Vengalattore, S. R. Leslie, J. Guzman, and D. M. 
Stamper-Kurn, Phys. Rev. Lett. 100, 170403 (2008). 
M. Vengalattore, J. Guzman, S. Leslie, F. Serwane, and 

D. M. Stamper-Kurn, arxiv:0901.3800. 

H. Saito, Y. Kawaguchi, and M. Ueda, Phys. Rev. A 76, 
043613 (2007). 

A. Lamacraft, Phys. Rev. Lett. 98, 160404 (2007). 



12 



[22] M. Uhlmann, R. Schiitzhold, and U. R. Fischer, Phys. 

Rev. Lett. 99, 120407 (2007). 
[23] B. Damski and W. H. Zurek, Phys. Rev. Lett. 99, 130402 

(2007). 

[24] S. Mukerjee, C. Xu, and J. E. Moore, Phys. Rev. B 76, 
104519 (2007). 

[25] G. I. Mias, N. R. Cooper, and S. M. Girvin, Phys. Rev. 

A 77, 023616 (2008). 
[26] R. W. Cherng, V. Gritsev, D. M. Stamper-Kurn, and 

E. Demler, Phys. Rev. Lett. 100, 180404 (2008). 
[27] T. Mizushima, K. Machida, and T. Kita, Phys. Rev. Lett. 

89, 030401 (2002). 
[28] E. J. Mueller, Phys. Rev. A 69, 033606 (2004). 
[29] T. Mizushima, N. Kobayashi, and K. Machida, Phys. 

Rev. A 70, 043613 (2004). 
[30] M. Stone, Phys. Rev. B 53, 16573 (1996). 
[31] C. V. Ciobanu, S. K. Yip, and T. L. Ho, Phys. Rev. A 

61, 033607 (2000). 
[32] R. Diener and T.-L. Ho, Phys. Rev. Lett. 96, 190405 



(2006). 

[33] L. Santos and T. Pfau, Phys. Rev. Lett. 96, 190404 
(2006). 

[34] E. Majorana, Nuovo Cimento 9, 43 (1932). 

[35] J. Schwinger, Quantum Theory of Angular Momentum 
(Academic, New York, 1994), edited by L. Biedenharn 
and H. Van Dam. 

[36] A. Auerbach, Interacting Electrons and Quantum Mag- 
netism (Springer, New York, 1994). 

[37] R. Barnett, A. Turner, and E. Demler, Phys. Rev. A 76, 
013605 (2007). 

[38] R. Rajaraman, Solitons and instantons (North Holland, 
1982). 

[39] R. Kamien, Rev. Mod. Phys. 74, 953 (2002). 
[40] U. A. Khawaja and H. Stoof, Nature 411, 918 (2001). 
[41] C. Bouchiat and G. W. Gibbons, J. Physique 49, 187 
(1988). 

[42] J. H. Hannay, J. Phys. A 31, 53 (1998). 



